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ABSTRACT 

The eigenvalue sensitivity for hydrodynamic stability operators is investigated. Classical 
matrix perturbation techniques as well as the concept of e-pseudoeigenvalues are applied to 
show that parts of the spectrum are highly sensitive to small perturbations. Applications are 
drawn from incompressible plane Couette, trailing line vortex flow and compressible Blasius 
boundary layer flow. Parametric studies indicate a monotonically increasing effect of the 
Reynolds number on the sensitivity. The phenomenon of eigenvalue sensitivity is due to the 
non-normality of the operators and their discrete matrix analogs and may be associated with 
large transient growth of the corresponding initial value problem. 
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1. Introduction 


Hydrodynamic stability theory has received a great deal of attention over the past 
decades (see e.g. [3] [13]) and still is a field of active and ongoing research. Its widespread 
applications range from basic fluid, dynamics to magnetohydrodynamics, from geo- 
physics to aeroelasticity, to name only a few. 

The starting point is a nonlinear evolution equation for a state vector describing the 
flow field. A basic state is assumed and a linearization of the nonlinear equation about 
this mean state results in a linear initial-boundary- value problem for the perturbations 
$ of the form d$/dt + i£$ = 0. The linear operator C depends in general on parameters 
describing the particular flow problem, such as Reynolds number, Mach number etc.. 
Linear stability is defined in an asymptotic manner for t — * oo. If for a given value of the 
flow parameters any (small) perturbation dies out, the system is said to be linearly stable 
for this particular parameter value. In the case of exponentially growing perturbations 
the system is linearly unstable and nonlinear effects will become significant once the 
initially small perturbations reach sufficient amplitudes. 

Eigenvalue analysis is the most dominant tool in determining the linear stability of a 
system. It is based on the assumption that the time dependence of the perturbations 
can be expressed in the form exp(-t'At). This reduces the linear initial-boundary-value 
problem to an eigenvalue problem A 'I* = £. 'I* . If the spectrum of £ consists of any 
eigenvalues whose imaginary part is positive, the system will support exponentially 
growing perturbations. If the spectrum of £ is confined to the lower complex half- 
plane, the system is linearly stable and small perturbations will decay for sufficiently 
large time. This technique of determining the linear stability has been used almost 
exclusively during the past decades. It is important, however, to stress that both linear 
stability and eigenvalue analysis are defined in an asymptotic setting ( t — > oo) and do 
not provide any information about the transient behavior of the operator £. 

In most cases, the operators involved in hydrodynamic stability problems are non- 
normal, i.e. although their set of eigenfunctions may be complete, the eigenfunctions do 
not form an orthogonal basis. Although this has been known for a long time [1], little 
attention has been paid to its consequences for the use of eigenvalues to describe the 
temporal behavior of infinitesimal disturbances. Non-normality may result in a large 
sensitivity of the eigenvalues to small perturbations of the operator. This has often been 
experienced as unsatisfactory convergence of parts of the spectrum and has been noted 
by several investigators [6] [14] [17] . As will be shown in this report, the sensitivity of the 
eigenvalues reveals more information about the short-time behavior of the operator than 
do the eigenvalues. Large transient growth is often associated with a high sensitivity 
of £ and can occur even if the eigenvalue spectrum predicts decay. This growth can 
take substantial values depending on the departure from normality of the operator or 
its discrete matrix analog. For this reason, the study of the sensitivity of hydrodynamic 
stability operators is as essential as the study of its eigenvalues. 

It has recently been shown [19] that parts of the Orr-Sommerfeld spectrum for plane 
Poiseuille flow are highly sensitive to small perturbations. Furthermore, even for Reynolds 
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numbers far less than the critical value of 5772, large transient growth of perturbation 
energy has been found. For a detailed analysis the reader is referred to [19]. 

In this paper, -we address the issue of eigenvalue sensitivity and discuss the use of eigen- 
value analysis for the investigation of transient non-mode-like behavior. Applications 
are drawn from different fields of fluid dynamics. Although limiting ourselves to exam- 
ples arising in fluid dynamics, we anticipate similar effects for linear stability problems 
that involve highly non-normal operators. 

This paper is organized as follows. Section 2 presents a brief review of some classical 
matrix perturbation concepts and introduces the idea of e-pseudoeigenvalues. Section 
3 applies the developed tools to the stability of both incompressible and compressible 
flows. Section 4 summarizes and concludes this paper. 

We would like to stress that this work is not concerned with roundoff error analysis. 
The focus of the present analysis is on the sensitivity of the spectrum as an effect of the 
non-normality of hydrodynamic stability operators. 


2. Basics of Matrix Perturbation Theory 

Matrix perturbation theory applied to eigenvalue problems is mainly concerned with the 
question of relating a small perturbation in the entries of a matrix A to the variation in 
its eigenvalues. Let E be a matrix with random entries and norm ||2£|| = £ <C 1 and let 
A be the perturbed matrix A + E. Given e and the (unperturbed) spectrum of A, how 
is the spectrum of A related to the spectrum of A ? One of the most straightforward 
estimates is given by linear perturbation theory [21]. 

Let A £ £ NxN b e a diagonalizable matrix, i.e. V~ l AV — diag(Ai , A 2 , ..., Ajv ). The 
Bauer-Fike theorem (see [5]) relates the maximum eigenvalue deviation (the spectral 
variation [21]) to the norm of the perturbation matrix in the following way. 

max min |Afc — Ay | < || V~ l EV\\ < k(V)\\E\1 (1) 

k j 

where A and A are eigenvalues of A and A, respectively, and k(V) = ||P -1 1| || W 1 1 denotes 
the condition number of the eigenvector matrix V. Throughout this study, we will 
use the 2-norm (Euclidean norm) denoted by ||a:|| = ||x ||2 = -\/(x, x ) with x € C N 
where (., .) stands for the associated inner product. The 2-norm generalizes to square 
matrices in an obvious way. Equation (1) states that, if each perturbed eigenvalue 
is associated with its closest unperturbed eigenvalue, the maximum deviation in the 
entire spectrum is bounded by the product of condition number of V and the norm of 
the perturbation matrix. The condition number n(V) = 1 if and only if A is normal, i.e. 
A h A — AA h where H denotes the conjugate transpose. Therefore, for normal matrices 
the eigenvalues of A will at most deviate from the eigenvalues of A by an amount of 
||£||. For non-normal matrices, k(V) can be substantially larger than 1 giving rise to 
0(1) deviations even for minute perturbations in the matrix entries. 
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The expression above bounds the largest deviation of an eigenvalue in the entire spec- 
trum as a function of ||£/||, but does not account for the sensitivity of individual eigenval- 
ues to small perturbations. A simple perturbation expansion for an individual eigenvalue 
yields [21] 

n H Ft, 

At = A. + + 0(||E|p) ||£|| « 1 (2) 

y"x k 

where x k and y k denote the right and left eigenvectors of A corresponding to the eigen- 
value A*. Neglecting higher order terms and taking the norm on both sides of equation 
(2), one obtains 


|A;t — Afc[ < s(A*)||-E|| 


( 3 ) 


with 


s(A fc ) 


MM 

\Vk*k\ 


( 4 ) 


Because of relation (3) the quantity s(X k ) is known as the condition of an eigenvalue [5] 
and can be interpreted as a measure of sensitivity to small perturbations in A. Again, 
for normal matrices s ( \ k ) = 1. 

Computing the condition s(Ajfc) of an eigenvalue A fc , is a straightforward procedure 
to obtain a first indication about parts of the spectrum that may experience a high 
sensitivity to small perturbations. 

In an effort to analyze the behavior of non-normal operators and matrices, Trefethen 
introduced the concept of e-pseudoeigenvalues [22] [23]. Although valid for operators 
and their discrete matrix analogs we will restrict ourselves to the matrix case only. 

Let A(A) be the spectrum of A € C NxN . The e-pseudospectrum A e (A) can be defined 
in two equivalent ways [22]. 


A e (A) = {zeC:\\(zI~Ar 1 \\>e- 1 } (5) 

or 

A e(A) = {z e C : z e A(A + E) for some E with ||I2|| < e}. (6) 

The first definition makes use of the norm of the resolvent R(z) = ( zl — A) -1 which 
is continuously defined in the complex plane with the exception of the spectrum A(A). 
As z approaches an eigenvalue of A, f?(z) will become singular and ||i?( 2 :)|| will tend to 
infinity. Therefore, the point spectrum of A will correspond to the locations z € C where 
the norm of the resolvent will become infinite. For finite but large resolvent norm, z 
will be defined as an e-pseudoeigenvalue of A with e = ||.R(z)|| -1 . 

The second definition makes use of the spectrum of a randomly perturbed matrix, z is 
considered an e-pseudoeigenvalue of A if it is an (exact) eigenvalue of a matrix A which 
is perturbed by a random matrix E of norm e. 
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Both definitions define the e-pseudospectrum as a region in the complex plane rather 
than a discrete set. Trefethen [22] [23] has shown that the e-pseudospectruin reveals more 
about the behavior of non-normal operators and matrices than does the spectrum alone. 
We illustrate this in the next section by calculating e-pseudoeigenvalues of highly non- 
normal matrices arising in hydrodynamic stability theory. For a proof of the equivalence 
of definitions (5) and (6) and for a more theoretical background on £-pseudoeigenvalues, 
the reader is referred to references [19] [22]. 

It is worth noting that the condition of an eigenvalue s(A^) can be used to bound the 
e-pseudospectrum of A. Combining (6) and (3), we find 


N 


A? £ |^J {X k -|- A s ( Afc ) e ] 


(7) 


where A r denotes a ball of radius r. This relation states that the e-pseudospectrum 
is enclosed in the union of disks around the eigenvalues A* with radius s(Ajt)e. It is 
important to keep in mind that relation (7) only holds for e sufficiently small such that 
higher order terms in equation (2) are negligible. Again, for normal matrices, that is 
■s(Ajt) = 1, we obtain a sharp bound on the e-pseudospectrum. 

An important implication of the sensitivity of the spectra is the possibility of large tran- 
sient growth of a disturbance governed by a non-normal operator. To demonstrate this, 
let us return to the linearized initial-boundary-value problem governing the evolution 
of infinitesimal perturbations, 

dd> 

+ iL<f> = 0 , ( 8 ) 

where <f> denotes the spatially discretized state vector and L now stands for the discretized 
matrix analog of the linear stability operator. Since the non-normality is an operator 
property, any consistent discretization technique for C will capture its behavior as long 
as sufficient resolution is ensured. The solution to (8) is given by 

4> = (f>o exp(-itL) (9) 

with (j > o as the initial condition <j)(t = 0). Applying eigenvalue analysis, it is straightfor- 
ward to show [23j that 

exp(-inf) < || exp(-i*.L)|| < «(F)exp(-fcrt), (10) 

where a stands for the largest imaginary part of the spectrum of L, i.e. the imaginary 
part of the least stable mode of L. 

If L is normal, then k(V) = 1 and the behavior of the initial-boundary- value problem 
is completely determined by the eigenvalues of L. In case L is highly non-normal, 
i.e. k(V) 1, the upper bound on || exp(— itL)\\ allows for transient growth, i.e. 
UW/WM li although the spectrum of L is confined to the lower complex half-plane. 
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For non-normal matrices, it is therefore necessary to compute the norm of the matrix 
exponential instead of relying on the spectrum of L to make predictions about the 
temporal transient behavior as described by the initial-boundary-value problem. 

A necessary and sufficient condition for the existence of transient growth can be derived 
in a straightforward manner from the evolution equation for the norm of the state vector. 
Assuming that the discretized linear initial value problem is given in the form of equation 
(8), the evolution equation for the norm of the state vector <f> reads 

JjIMI = (%•>’) + (*, %) = (-•'£*. « + W, = a»IW, «]. (a) 

where the relation ( f,g ) = (</,/)* has been utilized. 

It is evident from equation (11) that temporal growth of <f) measured in the Euclidean 
norm can occur if and only if the imaginary part of (L<f), <j>) is positive. The quantity 
( L(f> , 4>) is known as the field of values or, equivalently, the numerical range of L and is 
defined as 


F{L) = {z : z = ( L4 >, 4>) with ||^|| = 1}. 

It is clear that the field of values of L contains the spectrum of L (see [8]). It can further 
be shown (e.g. [22]) that for normal matrices, E(L) is the convex hull of the spectrum 
whereas for non-normal matrices the field of values can be a considerably larger set 
than the spectrum. Especially cases in which the field of values E(L) reaches into the 
unstable domain, although the spectrum of L is confined to the stable half-plane, are of 
particular interest. For these cases, transient growth followed by asymptotic decay can 
be predicted. 

The boundary of the field of values gives an immediate indication for transient growth 
and can easily be calculated using standard techniques [8]. It should be kept in mind, 
however, that the field of values JF(L) is depends on the choice of the scalar product 
and is an indicator for growth measured in its associated norm. 

3. Application to Hydrodynamic Stability Problems 

In this section we apply the tools developed in the previous chapter to study the sen- 
sitivity of an assortment of matrices resulting from discretized hydrodynamic stability 
operators. We thereby use the following qualitative and quantitative techniques to 
illustrate the properties of the spectra under consideration. 

• Following the second definition of £-pseudospectra we display a superposition of 
100 spectra of randomly perturbed matrices A = A + E. The complex entries of 
the matrix E are drawn randomly from a normal distribution with mean zero and 
standard deviation one. The norm of the matrix E is chosen to be e. The plot then 
shows £-pseudoeigenvalues in a “statistical” fashion and gives a first impression of 
the sensitive parts of the spectrum under consideration. 
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• For a few eigenvalues the condition defined by (4) is tabulated which gives an upper 
bound for the eigenvalue deviation (equation (3)) for sufficiently small perturba- 
tions. Although the bound might not be sharp, it gives insight into the relative 
sensitivity of eigenvalues to small perturbations. 

• An illustrative means to understand the spectrum of a non-normal operator is a 
contour plot of constant resolvent norm || (zI—A)~ l || = 1/e for various e. It displays 
closed, nested curves that confine the eigenvalues of the matrix A if perturbed by 
a random matrix of norm e. The numerical realization of this kind of plots is 
rather time consuming and involves a singular value decomposition on grid points 
discretizing the region of interest in the complex plane (see [19] [22]). 

We consider three different flow types to demonstrate spectral properties of non- 
normal operators. The first example is incompressible, plane Couette flow, i.e. the flow 
between infinite plates moving in opposite direction. Another application comes from 
trailing line vortex flow. Both inviscid and viscous disturbances are considered. The 
final example is taken from compressible Blasius boundary layer flow, i.e. flow in a 
semi-infinite domain. 

We would like to emphasize that, as long as sufficient resolution is ensured, the 
method of discretizing the stability operator does not influence the qualitative results 
of this study. The sensitivity of the spectrum is a property of the operator that will be 
captured by its discrete matrix analog. 

3.1. Plane Couette Flow 

Plane, incompressible Couette flow, i.e. flow between infinite parallel plates moving 
in opposite direction, can be considered as one of the canonical flow situations. The 
linear mean velocity profile greatly simplifies theoretical investigations, which made this 
type of flow particularly accessible to analytical methods. Herron [7] showed that plane 
Couette flow is linearly stable for all Reynolds numbers. 

The problem of transient growth in plane Couette flow has been studied previously by 
Farrell [4] using variational principles and recently by Reddy & Henningson [18] making 
use of ^-pseudospectra. In this presentation, it is included for completeness and is meant 
as an introductory example to familiarize the reader with the analysis of non-normal 
matrices and their associated characteristics. 

The equation governing the evolution of infinitesimal, two-dimensional disturbances 
in plane Couette flow of the form v(x,y,t) = t)(y, t) exp(iax) can be written as 

(12a) 

v(±l) = Dv(± 1) = 0, (12 b) 

with 

C = (D 2 - a 2 ) -1 [aU(D 2 - a 2 ) - ^-(D 2 - a 2 ) 2 ] (12 c) 

iR 


dv 

dt 


+ iCv = 0 
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where v denotes the Fourier transformed normal velocity, a represents the streamwise 
wavenumber, U = y stands for the mean velocity profile and D denotes differentiation 
in the wall normal direction y. All flow quantities have been non-dimensionalized by 
the velocity of the moving plates Uo and the channel half-height h and the Reynolds 
number R is based on these scaling parameters, i.e. R = Uoh/v. The equation above 
is related to the Orr-Sommerfelcl equation which has received a great deal of attention 
over the past decades [2]. 

The operator C is spatially discretized using a hybrid Chebyshev pseudospectral tech- 
nique and the discretized operator is denoted by L. 
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Fig.l. Pseudospectra for plane Couette How for a = 1 and R = 1500. (a) Superpo- 
sition of 100 spectra perturbed by random matrices of norm 10 -6 . The unperturbed 
spectrum is displayed by the larger square symbols, (b) Contour plot of constant resol- 
vent norm. The contours represent levels from 10 -8 ( innermost contour) to 10 -2 (outer 
contour). 
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Figure la shows the superposition of 100 spectra of L for a = 1 and R = 1500 per- 
turbed by random matrices of norm 10 — 6 . The unperturbed spectrum is shown by larger 
squares. It is observed that a perturbation of norm 10“ 6 can cause some eigenvalues 
to move by a distance of order one. This sensitivity is highest for eigenvalues at the 
intersection point of the three eigenvalue branches. This observation is also apparent 
from Table 1 which lists selected eigenvalues together with their condition as defined by 
(4). It can be seen that the condition is largest for the eigenvalue in the junction point 
whereas the least stable mode has a rather low condition, which is also reflected in the 
low sensitivity of the least stable modes in Figure la. Figure lb, which displays the 
resolvent norm contours together with the unperturbed spectrum, gives the complete 
picture of the sensitivity properties of the matrix L for the specified parameter values. 
It is seen that even a perturbation of 10 -8 has an appreciable effect on the location of 
the eigenvalues near the intersection point of the three branches. In order to move the 
least stable modes a distance of order one, a perturbation of norm 10 2 is needed which 
renders them rather insensitive. It should be noted that the 10 2 -pseudospectrum ex- 
tends into the upper half-plane which has implications on the possibility of transient 
growth. For more information on this issue the reader is referred to [19]. 


eigenvalue 

condition 

A r 

A, 


0.723771 

-0.259583 

0.178337 10 2 

0.653132 

-0.102713 

0.289772 10 2 

0.336625 

-0.457634 

0.328836 10 5 

-0.516156 10~ 5 

-1.42610 

0.833198 10 6 

0.218617 10“ 3 

-0.702198 

0.167807 10 9 


Table. 1. Selected eigenvalues and conditions for plane Couette flow for a — 1 and 
R = 1500. 

It is interesting to probe the transient temporal behavior of L. To this end, we 
calculate the norm of the matrix exponential of L for a = 1 and R = 1500. Figure 2 
shows || exp(-itL)\\ versus time. Large transient growth up to about ||^||/||^o|| ~ 20 
is observed before exponential decay predicted by the eigenvalues sets in. The reason 
for this growth lies in the non-orthogonality of the set of eigenvectors of L, which in 
turn is an effect of the non-normality of the respective stability operator. Figure 2 is 
a striking and instructive demonstration that predictions about the temporal behavior 
based on the eigenvalues of L can be misleading and that the non-mode-like behavior 
of the stability operator is not captured by its spectrum. 

As explained in section 2, the field of values of the linear stability operator for plane 
Couette flow can be used as an indicator for transient growth. For the parameter values 
chosen above, the maximum imaginary part of the field of values has been calculated 
to be 0.557 which is a measure for the largest growth rate possible. 
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Although measuring the growth in the 2-norm does reveal the transient behavior of 
the operator under investigation, a weighted norm has to be used to draw conclusions 
about the amplification of physical flow quantities such as perturbation energy; for 
details the reader is referred to [19]. 



Fig.2. Plot of || exp(— HL)\\ versus time for plane Couette flow at a = 1 and R = 
1500. 

The dependence of the sensitivity on the Reynolds number is shown in Table 2. For 
this case, the sensitivity has been defined as — logioGK-i I ~ £) -, |l)> he. the negative 
exponent of the resolvent norm evaluated at selected locations z } in the complex plane. 
Table 2 shows a monotonic increase of the sensitivity with increasing Reynolds number 
which is more pronounced for locations near the intersection point of the three eigenvalue 
branches than near the least stable modes. 


Reynolds number 

z = 0.7- i 0.2 

sensitivity 
z = 0.4 - i 0.4 

z = - i 0.4 

500 

1.772 

3.092 

4.544 

1000 

2.095 

4.343 

5.277 

1500 

2.752 

4.975 

6.770 

2000 

2.850 

5.790 

7.232 

2500 

3.020 

6.434 

8.446 


Table. 2. Dependence of the sensitivity on the Reynolds number at selected locations in 
the complex plane. 
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3.2. Trailing Line Vortex Flow 


Stability analysis of swirling flows is an active area of research. One particular mean 
flow which has received much attention in recent years is the trailing line (Batche- 
lor) vortex model([9]-[12],[l6]). In cylindrical coordinates (r,9,z) the base flow for the 
columnar vortex is given by 


<7 = 0, V = -(1 — e~ r ’ 2 ), W = W Q + c- r \ (13) 

r 

where <7, V , and W are the radial, tangential, and axial velocities, respectively. The 
swirl parameter q is related to the ratio of maximum azimuthal velocity to maximum 
axial velocity and Wo represents the undisturbed outer axial flow. Without loss of 
generality, throughout this section we assume Wo = 1 . 

The above vortical mean flow is perturbed by three-dimensional disturbances of the 

type 


[u, v, w, p] = [ iu(r , t), v(r, t ), w(r, t), p(r, *)] e l ( az + 11 ^ , (14) 

where a is the axial wavenumber and n is the azimuthal wavenumber which takes on 
integer values only. As in the previous section, the linearized governing equations (see 
[ 9 ], [ 10 ]) are arranged in the form d(j>/dt + i£(j> — 0 , where 4> is the state vector 


<f> = [u(r, t ), v(r, t), w(r , t), p(r, t)] . 
The operator C is given by 


C = 


( Jn 

,_. Tr 2 in V 

(m ~w + 7 ) 

DW 

-i(D+b 

7 r 


2V 2 in 
r Rr 2 


I 22 

0 


m 
7 r 


) 0 

0 

hi 
ia 

7 


-d\ 


n 

r 

a 


— — 0 


with 


iD 


nV 


i , n 2 + 1 


ln = h2 = ^(D + l) + ^ + aW-^^+^) 


R 


(15) 


(16a) 


and 

'»» = ,n + 
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In equation ( 16 a), D represents the radial derivative and R is the Reynolds number 
based on the viscous core radius and the axial velocity excess (or deficit). 

The boundary conditions far from the centerline of the vortex are 

u( oo) = 6(oo) = ui(oo) = 0. ( 166 ) 

For single- valued and smooth solutions, the boundary conditions on the centerline of 
the vortex are 

if n = 0, 6(0) = 6(0) = 0, w( 0) and p(0) must be finite, 

if n = ±1, u(0) ± 6(0) = 0, D6(0) = 0, u>(0) = p(0) = 0 ( 16 c) 

if |n| > 1, 6(0) = 6(0) = i6(0) = p( 0) = 0. 

The form of the continuity equation poses a problem as it does not contain a time 
dependent term and therefore cannot be cast into an eigenvalue problem. Although row 
and column operations can be used to solve this issue, a more efficient method, used 
in equation ( 16 ), is to artificially introduce the term 'ydp/dt to the continuity equation 
where the constant 7 can be thought of as a weight parameter. This technique is known 
as the artificial compressibility method. For swirling flows, Khorrami, Malik & Ash[ 9 ] 
showed that for very small values of the constant 7, the physical eigenvalues remain 
unaffected and the spectrum is computed efficiently. 

The operator C is spatially discretized using a Chebyshev spectral collocation tech- 
nique. The details of the discretization are given in [ 9 ] and are omitted here. 

Two distinct stability characteristics have resulted in a classification of swirling sta- 
bility modes into inviscid and viscous perturbations. Modes with negative azimuthal 
wavenumber are referred to as inviscid modes because their stability characteristics are 
determined by inviscid phenomena. The so-called unstable viscous modes are dominated 
by diffusive effects and include eigensolutions for azimuthal wavenumbers of n = 0,1. 
For large Reynolds numbers R >• 1, the growth rate of the viscous disturbances is 
inversely proportional to the Reynolds number, i.e. in the limit of infinite Reynolds 
number the viscous modes become neutrally stable. Based on previous investigations 
[ 10 ] [ 1 1 ] , the trailing line vortex represented by equation ( 13 ) is unstable with respect to 
both inviscid and viscous disturbances. 

Due to space limitations, in this study we restrict our attention to a single azimuthal 
Fourier mode taken from the inviscid and viscous family. We chose parameter values 
that have been found to be the most interesting/challenging from a computational point 
of view. 

Figure 3 displays the pseudospectrum for the inviscid case of trailing line vortex 
flow. The parameters are chosen as follows: a = 0 . 4 , n = — 1 ,R = lO 4 ,^ = 0.4 and 
7 = 10 -8 . Figure 3 a shows the e-pseudospectrum for e = 10 -5 . Sensitive eigenvalues 
are observed close to the neutral line whereas the least stable modes show little effect 
to the perturbations imposed. Figure 3b displays the resolvent contours for the same 
case where the contour levels range from 10 ~ 6 to 10 -4 . 
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real 


Fig. 3. Pseudospectra for trailing line vortex Bow ( inviscid case) for cv 0.4, n 
— = 10 A ,q — 0.4 and 7 = 10 -8 . (a) Superposition of 100 spectra perturbed by 

random matrices of norm 10~ 5 . The unperturbed spectrum is displayed by the larger 
square symbols, (b) Contour plot of constant resolvent norm . The contours represent 
levels from 10 -6 (innermost contour) to 10~ 4 (outer contour). Note the different scaling 
for the real and imaginary axis. 
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Table 3 lists the condition of selected eigenvalues and confirms the relative sensitivity 
shown in Figure 3. 

As was the case for Couette flow, the sensitivity of the eigenvalues increases mono- 
tonically. Table 4 presents the sensitivity evaluated at three different locations in the 
complex plane as a function of Reynolds number. In all three cases we experience an 
increase in sensitivity with increasing Reynolds number. 


€ 

A r 

iigenvalue 

A, 

condition 

0.436626 

0.383041 

0.366548 

0.307724 

0.356920 

0.182285 
0.914185 10 _1 
-0.735235 10" 1 
0.120135 10 _1 
0.248804 IQ" 1 

0.367540 10 3 
0.669828 10 3 
0.602626 10 3 
0.830603 10 2 
0.400462 10 4 


Table. 3. Selected eigenvalues and conditions for swirling How (inviscid case). Parameters 
as in caption of Figure 3. 


Reynolds number 

z = 0.42 + i 0.12 

sensitivity 
z = 0.375 - i 0.03 

z = 0.365 + i 0.02 

1000 

4.076 

4.176 

5.236 

2000 

4.092 

4.319 

5.279 

5000 

4.111 

4.683 

5.322 

10000 

4.121 

5.029 

5.635 

20000 

4.128 

5.228 

5.968 


Table. 4. Dependence of the sensitivity on the Reynolds number ( inviscid case). 


The pseudospectrum for the viscous case of the trailing line vortex flow is shown 
in Figure 4 for a = 0.3, n = 1 ,R= 10 4 ,<? = 0.4 and 7 = 10 -9 . Figure 4 a reveals 
that a perturbation of norm 10 -6 results in a deviation of the eigenvalues orders of 
magnitude larger than the perturbation imposed. Figure 4 b displays the contours of 
constant resolvent norm with levels ranging from 10 -8 to 10 -5 . It is interesting to note 
that the unstable eigenvalue lies within the 10 _ 8 -contour. 
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real 

Fie:. 4. Pseudospectra for trailing line vortex flow ( viscous case) for a = 0.3, n = 
1 r — io 4 o = 04 and 7 = 10 -9 . (a) Superposition of 100 spectra perturbed by 

random matrices of norm 1 (T 6 . The unperturbed spectrum is displayed by the larger 
square symbols, (b) Contour plot of constant resolvent norm. The contours represent 
levels from 1(T 8 (innermost contour) to 1CT 5 ( outer contour). Note the different scaling 
for the real and imaginary axis. 
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Table 5, listing the condition for assorted eigenvalues, emphasizes this fact as the 
unstable mode shows the largest eigenvalue condition. A parameter study, summarized 
in Table 6, shows that the sensitivity of the spectrum, especially near the unstable 
mode, grows quite appreciably. 


€ 

A y 

igcn value 

A t 

condition 

0.380991 

0.350173 

0.323453 

0.322607 

0.310078 

-0.527545 10~ 2 
-0.339581 10~ 2 
0.918111 10~ 4 
-0.155546 10~ 2 
-0.655544 10“ 3 

0.231711 10 3 
0.522167 10 3 
0.159909 10 6 
0.150151 10 6 
0.213876 10 4 


Table. 5, Selected eigenvalues and conditions for swirling Bow (viscous case). Parameters 
as in caption of Figure 4. 


Reynolds number 

z = 0.38 - i 0.004 

sensitivity 
z = 0.35 - i 0.002 

z = 0.322 - i 0.001 

30 

4.089 

4.391 

4.970 

60 

4.381 

4.755 

5.490 

200 

4.864 

5.367 

6.469 

500 

5.163 

5.728 

7.087 

1000 

5.362 

5.946 

7.408 

3000 

5.657 

6.251 

7.777 

10000 

5.772 

6.407 

8.220 


Table. 6. Dependence of the sensitivity on the Reynolds number 


Sensitive eigenvalues close to the neutral line can result in large transient growth 
of the perturbation amplitude. In order to show that, we calculated the norm of the 
matrix exponential as described in section 2. Figure 7 shows the amplification factor 
UW/UoW versus time. A large transient growth is observed which only for large times 
is governed by the unstable eigenvalue. 

The maximum imaginary part of the field of values has been determined to be 24.78 
which confirms the large growth of the state vector norm as plotted in Fig. 5. 
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Fig. 5. Plot of || exp(-ziL)|| versus time for trailing line vortex flow (viscous case) for 
parameters given in caption of Figure 4. 


3.3. Compressible Blasius Boundary Layer Flow 

In this section we investigate the spatial evolution of infinitesimal disturbances as 
opposed to the temporal evolution in the previous examples. 

Consider compressible boundary-layer flow past a flat-plate. It is convenient to write 
the basic state using the similarity transformation 

* - M (plp ' )iy (17) 

where ,r, y are the coordinates parallel and normal to the plate, p e is the boundary- 
layer edge density, u e the streamwise velocity and p e the dynamic viscosity. The gov- 
erning equations for the mean flow can be written in a similarity form as 

(cry + if " = ° ( l8 ) 

(dig 1 + f")' + fg 1 = Q (^) 


where 


f' = u/u e , C = pp/p e p e 
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g = H/H„ a, = c/< 7 , a 2 = — - (1 - -)c 

and A/ is the Mach number, 7 the ratio of specific heats, a the Prandtl number, and 
H the total enthalpy. 

The compressible boundary- layer flow is assumed to be locally parallel and is per- 
turbed by a harmonic disturbance of the form 

$ = $(x, y)e l (P z ~ (20) 

where (3 is the spanwise wavenumber and cj is the disturbance frequency. We consider 
the spatial stability in which disturbances grow in the streamwise direction x. The 
spatial eigenvalue problem is nonlinear due to the streamwise viscous diffusion ($ xr ) 
term. It was argued in [15] that the eigenvalue problem can be linearized by dropping 
the <f> xx term since the instability in supersonic flows is predominantly inviscid. Later 
work suggests that this approximation is valid even at lower Mach numbers and is also 
consistent with the triple-deck description of the Tollmien-Schlichting wave structure. 
Hence, the governing compressible linear stability equations can be written as 


3 $ 

— + - 0 

ox 


( 21 ) 


where <5 is a five-element vector defined 
and the linear operator C is of the form 


by [$1, $2, 3*3, $4, $5} = , 


C — A — + B — — t C (22) 

dy 2 dy 

where A , B and C are 5x5 matrices whose elements are related to the matrices given 
in [15]. 

Equation ( 21 ) is solved subject to the boundary conditions 

f li = 4*2 = d ) 4 = $5=0 at y = 0, (23a) 

4>i , $ 2 , $ 4 , $5 — * ► 0 as y — > oo. (236) 

In the spatial analysis, the x-dependence of the disturbance is expressed in the form 
exp(iAx) which results in the eigenvalue problem — A$ = £$. Disturbances are expected 
to grow in x if Im(A) < 0. 

The operator C is discretized by a Chebyshev spectral collocation technique described 
in [15]. 

The pseudospectrum for compressible boundary layer flow is computed for a Reynolds 
number R = 1000, where R = \Jp^u^xjy e . In Fig. 6 , it is displayed for the parameter 
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combination M = 2,/3 = 0.08 and u> — 0.02. Figure 6a shows a superposition of 100 
spectra perturbed by random matrices of norm 10 -3 . Although the sensitivity of the 
eigenvalues is less pronounced than in the previous cases, parts of the spectrum do show 
eigenvalues that deviate orders of magnitude more than the perturbation that caused the 
deviation. Especially the least stable mode and the continuous entropy branch exhibit 
a marked response to the perturbations as displayed in Figure 6b. This observation is 
again confirmed by calculating the condition of the respective eigenvalues (Table 7). 


eig< 

A r 

snvalue 

A t 

condition 

0.407748 KT 1 
0.200015 10' 1 
0.200007 10 _1 
0.529917 10" 1 
0.108886 

-0.239187 10~ 2 
0.115666 10- 1 
0.460632 10- 1 
0.338088 lO" 1 
0.178297 10~ 3 

0.210794 10 2 
0.204110 10 1 
0.190301 10 1 
0.425356 10 2 
0.228352 10 1 


Table. 7. Selected eigenvalues and conditions for compressible boundary layer flow. Pa- 
rameters as in caption of Figure 6. 
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Fig. 6. Pseudospectra for compressible boundary layer How for M = 2 . fi = 0 . 08 , R — 
1000 and u> = 0 . 02 . (a) Superposition of 100 spectra perturbed by random matrices of 
norm 10 -3 . (b) Contour plot of constant resolvent norm. The contours represent levels 
from 10 -4 (innermost contour) to 10~ 2 (outer contour). The spectrum is displayed by 
the square symbols. Note the different scaling for the real and imaginary axis. 
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4. Summary and Conclusions 


A sensitivity study of hydrodynamic stability operators has been conducted to in- 
vestigate the response of eigenvalues to small perturbations. Classical techniques from 
matrix perturbation theory as well as the concept of e-pseudoeigenvalues has been ap- 
plied to obtain qualitative and quantitative results. Applications involved incompress- 
ible plane Couette flow, incompressible trailing line vortex flow and compressible Blasius 
boundary layer flow. All examples showed sensitive eigenvalues that deviated from their 
unperturbed location by orders of magnitude more than the norm of the perturbations 
imposed on the matrix entries. 

In the case of Couette flow, the most sensitive eigenvalues have been found in the 
junction of the three eigenvalue branches whereas the least stable eigenvalues showed a 
rather small response to small random perturbations. The sensitivity of the spectrum 
measured by the norm of the resolvent evaluated at specific points in the complex plane 
increased monotonically with Reynolds number. 

For the case of trailing line vortex flow, viscous and inviscid perturbations have to be 
distinguished. For inviscid perturbations, again the unstable mode seems to be rather 
insensitive and the regions of high sensitivity are located close to the real axis. The 
case of viscous disturbances imposed on the vortex flow, however, showed a dramatic 
effect of small perturbations on the location of the eigenvalues. The only unstable 
mode experienced the highest sensitivity implying the loss of its predictive character 
because of the dependence on very small perturbations. As in Couette flow, in both 
cases of trailing line vortex flow the sensitivity increased monotonically with Reynolds 
number. Investigations of the norm of the matrix exponential revealed large transient 
growth which resulted in a growth rate superior to the one predicted by the unstable 
eigenvalue. The field of values of the linear operator gave a necessary and sufficient 
condition for transient growth as well as a measure for the maximum possible growth 
rate. 

Although the compressible Blasius boundary layer spectrum appeared to be less sensi- 
tive to small perturbations than the previous cases, only a small perturbation is needed 
to find eigenvalues stemming from the continuous entropy branch in the unstable half- 
plane. It is conjectured that the sensitivity will increase for higher Reynolds and/or 
Mach numbers. 

One of the main interest of hydrodynamic stability theory is the understanding of 
the temporal behavior of small perturbations superimposed on a steady flow field. This 
behavior is governed by the initial-boundary-value problem stated in section 1 of this 
presentation. The analysis can, therefore, be restricted to understanding the behavior 
of the operator C. 

There are two distinct regimes to be considered: the asymptotic regime of large times 
(f — > oo) and the transient regime. In the first regime, it has been shown that the 
eigenvalues of C are sufficient to completely describe the behavior of C. In the second 
regime, however, eigenvalues do not suffice to describe the temporal behavior since it 
may not be mode-like. Hence, there is a need to replace eigenvalue analysis by a more 
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effective tool to analyze operators which can experience large growth due to their non- 
normality during this transient period, e-pseudoeigenvalues, i.e. eigenvalues of the 
slightly perturbed operator, can be used to extract informations about the operator in 
its transient phase and to “visualize its behavior”. 

As the transient growth can be quite substantial, the nonlinear regime may be reached 
during the transient phase, rather than by a mechanism based on exponential growth. 
As has been shown [20], the growth mechanism stemming from the non-normality of the 
linear operator plays a dominant role during the early stages of transition to turbulence 
in parallel shear flow. For incompressible flows it is only this linear mechanism that 
can provide an increase in total perturbation energy. The implication of growth due 
to non-normality is therefore of primary importance in hydrodynamic stability theory. 
A necessary condition for growth is the non-normality of the stability operator and 
sufficient conditions can be derived (together with bounds on the maximum growth 
possible) from the e-pseudospectrum [19]. 

For these reasons, the study of the e-pseudospectrum is essential when encountering 
highly non-normal operators or matrices. 
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